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A stochastic magnetic field in the early Universe will produce anisotropics in the temperature and 
polarization of the cosmic microwave background. We derive analytic expressions for the microwave 
background temperature and polarization power spectra induced by vector and tensor perturbations 
from a power-law magnetic field. For a scale-invariant stochastic magnetic field smoothed over a 
comoving scale of 1 Mpc, the MAP satellite has the potential to constrain the comoving mean-field 
amplitude to be no greater than approximately 2 x 10~® G. Limits improve as the power-law slope 
increases: for causally-generated power-law magnetic fields, the comoving mean-field amplitude has 
an upper bound of approximately 4 x 10~^^ G. Such constraints will surpass all current limits on 
galactic-scale primordial stochastic magnetic fields at decoupling. 
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I. INTRODUCTION 

Magnetic fields of strength are ubiquitous in galaxies 0] and clusters of galaxies [0. The origin of these fields, 
however, remains an outstanding problem in cosmology. It is usually postulated that these /iG fields grew either 
via some magnetohydrodynamical (MHD) dynamo mechanism or via adiabatic compression of a primordial 
magnetic field during the collapse of a protogalactic cloud . A MHD dynamo requires tiny seed magnetic fields of 
comoving amplitude 10"'^° G in conventional GDM-like cosmological models or even as tiny as 10^^° G in a Universe 
with a non-zero cosmological constant Has suggested by recent measurements of type la snpernovae and the 

microwave background power spectrum TlU-lIa] . On the other hand, the adiabatic compression scenario requires a far 
larger primordial seed field with comoving amplitude of 10^^ G to 10^^" G. 

Persistent questions about the effectiveness of MHD dynamos ]l6|-p3t together with the observation of /iG magnetic 
fields in high-redshift galaxies raise the possibility of a significant primordial magnetic field in galaxies and clusters 
of galaxies. The origin of such a magnetic field remains a mystery. Essentially all viable magnetogenesis mechanisms 
incorporate speculative ideas in high-energy theory, including (among others) inflation p^-p9|, electroweak ^0 31 



or QCD 1^,^ phase transitions, charge asymmetry ||3j], or a ferromagnetic Yang-Mills vacuum state As the 
properties of the primordial magnetic field predicted varies among these mechanisms, future detections of a primordial 
magnetic field may aid us in identifying the correct magnetogenesis mechanism. On the cosmological front, a primordial 
magnetic field may have affected early- Universe processes such as phase transitions, baryogenesis, and nucleosynthesis 
(see for a review). Relic magnetic fields could provide a direct source of information about these processes. A 
primordial ma gne tic field may also have infiuenced structure formation via contributing to density perturbations on 
galactic scales |37H40|] and preserving magnetic energy in Alfven modes on scales below the Silk damping scale during 
recombination | |4l| , ]42[ . In short, significant primordial magnetic fields would impact both cosmology and particle 
physics. 

The presence of a magnetic field in the early Universe affects the evolution of metric perturbations, and as a result, 
produces temperature and polarization anisotropies in the cosmic microwave background (CMB). High-resolution 
measurements of the microwave background provide a clean and model-independent test for primordial magnetic 
fields. We demonstrate in this paper that fields large enough to result in observed fields via adiabatic compression will 
likely leave observable and distinctive fluctuations in the various power spectra of microwave background temperature 
and polarization fluctuations. 

Substantial progress has been made in understanding the effects of a primordial magnetic field on the GMB. The 
Far Infrared Absolute Spectrophotometer (FIRAS) upper limits on chemical potential /i and Compton y distortions 
in the CMB blackbody constrain the present strength of the magnetic field with comoving coherence length between 
400 pc and 0.6 Mpc to be < 3 x 10~^ G |^]. The case of a homogeneous magnetic field has been considered 
by several authors. The best current constraint on the primordial homogeneous magnetic field strength is Bq < 
3.4 X 10~^(Oo'i5o)^^^ G (/i5o is the present Hubble constant in units of 50 km s~^ Mpc~^), obtained by doing 
statistical analysis on the 4-year Cosmic Background Explorer (COBE) data for temperature patterns of a Bianchi 
type VII anisotropic spacetime [H. A primordial homogeneous magnetic field can produce distortions of the CMB 
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acoustic peaks via fast magnetosonic waves [ p5[ ; meanwhile, Alfven wave excitations can amplify vector perturbations 
and induce additional correlations in temperature multipole moments p6|] . It is shown in Ref . that a primordial 
homogeneous magnetic field of present strength 10~^ G at decoupling can induce a measurable Faraday rotation in 
the CMB polarization of 1° at a frequency of 30 GHz. Additional CMB polarization effects arising from a primordial 
homogeneous magnetic field via Faraday rotation include a parity-odd cross correlation between temperature and 
polarization anisotropics and the depolarization of the original CMB polarization |^ , which leads to a reduction 
in the damping of temperature anisotropics on small angular scales. 

The case of a stochastic magnetic field is perhaps more realistic, because such fields are observed within galaxy 
clusters ||50|-|53| and predicted by all causal magnetogenesis mechanisms Q. Some numerical estimates of CMB 
temperature and polarization power spectra from density perturbations induced by a primordial stochastic magnetic 
field are presented in Ref. p^ , whereas corresponding analytic estimates, though somewhat crude and valid only 
for temperature anisotropics on large angular scales, are given in Ref. |5^. Effects of Alfven waves induced by a 
primordial stochastic magnetic field on CMB temperature and B-type polarization anisotropics are considered in 
Refs. 56 1 and jS^ respectively. Finally, a primordial stochastic magnetic field also generates gravitational waves; the 
resulting tensor CMB temperature power spectrum is given in Ref. [Esl . 

Although a variety of effects of a primordial stochastic magnetic field on the CMB have been investigated, the 
results are fragmented and a systematic approach is lacking. Besides the temperature power spectrum from tensor 
perturbations given in Ref. , no other CMB power spectra have been derived. We consider a statistically homoge- 
neous and isotropic stochastic magnetic field with a power-law power spectrum, generated at some early epoch of the 
radiation-dominated Universe. Based on the computational techniques in Ref. [ p8[ and the total angular momentum 
method for calculating CMB anisotropics introduced by Hu and White ||5^ , we have completed a comprehensive and 
unified analytic calculation of all types of CMB power spectra arising from a primordial stochastic magnetic field. 
This paper focuses on the induced vector and tensor perturbations. A primordial magnetic field acts as a continuous 
source of vorticity until decoupling and gravitational radiation until matter-radiation equality. The resulting vector 
and tensor perturbations are one of the few cosmological sources of B-type polarization [^,0 , along with primordial 
tensor perturbations and gravitational lensing of the CMB Scalar perturbations induce CMB aniso tropics 
smaller than those from vector and tensor perturbations and hence will not be considered here (see Sec. VIII). 

In Sec. H we derive the power spectrum for a primordial stochastic magnetic field. We then project the vector 
and tensor pieces from the electromagnetic stress-energy tensor, from which we calculate their two-point correlation 
functions and derive their isotropic spectra. Details of the derivation of the vector isotropic spectrum are presented in 
the Appendix. Section III computes the magnetic damping scales of the induced vector perturbations at decoupling 



and tensor perturbations at matter-radiation equality. In Sec. IV we review the vector and tensor contributions 
to the metric tensor and give their corresponding evolution equations. We obtain solutions to these equations, 
which can be expressed as functions of the magnetic-induced isotropic spectra derived in Sec. |l[ Using the total 
angular momentum method of Ref. we compute analytically the CMB pow er sp ectra for temperature in Sec. 
polarization in Sec. VI, and the temperature-polarization cross correlation in Sec. VII. Section VIII| concludes with the 
physical interpretation of these results and a discussion of current and future limits on primordial magnetic fields from 
the microwave background. For the vector perturbations, the B-type is slightly larger than the E-type polarization 
power spectrum, whereas the E-type polarization and the cross-correlation power spectra are approximately identical. 
For the tensor perturbations, the polarization power spectra are actually comparable to the temperature power 



spectrum for n > —3/2, where n is the magnetic field power-law spectral index. As we will show in Sec. HI, the 
tensor perturbations are damped on smaller scales than the vector perturbations. The magnetic cutoff wavenumber 
determines the overall amplitude of the CMB power spectra, so the tensor-induced CMB anisotropics will be larger 
than the vector anisotropics for n > —3/2. For a scale-invariant stochastic magnetic field smoothed over a comoving 
scale of 1 Mpc, near-future microwave background temperature measurements will constrain the comoving mean-field 
amplitude to be no greater than approximately 2x 10~^ G. Limits improve as n increases: for causally-generated power- 
law magnetic fields with n > 2, the comoving mean-field amplitude will soon have an upper bound of approximately 
4 X 10^^'^ G. These will be the strongest current constraints on galactic-scale primordial stochastic magnetic fields at 
decoupling. Eventually, precision measurements of the microwave background temperature and polarization will give 
significantly stronger constraints. 

In this paper, we focus on the induced CMB anisotropics for I < 500, where the analysis is relatively clean, simple, 
and free from complications arising from the last-scattering microphysics. For simplicity, we consider the case of a 
flat Universe with a vanishing cosmological constant. We employ the following notational conventions: a is the scale 
factor, 77 is the conformal time, overdots are derivatives with respect to 77, and subscripts denote the present time. 
We set c = 1 and normalize the scale factor to unity today. As usual, Greek indices run from to 3 and Latin ones 
from 1 to 3. All calculations are done in Fourier space, unless real-space dependence is indicated explicitly (as in 
Sec. All magnetic field amplitudes are comoving values, unless an explicit time dependence is displayed. 
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II. MAGNETIC POWER SPECTRUM AND CORRELATION FUNCTIONS 



Consider a primordial stochastic magnetic field created at some specific moment during the radiation-dominated 
epoch. The energy density of the magnetic field is treated as a first-order perturbation to a flat Friedmann-Robertson- 
Walker (FRW) background cosmology. In other words, we do not decompose the magnetic field into a large homoge- 
neous component and a small fluctuating piece as in most cases in the literature. Within the linear approximation, the 
magnetic field evolves as a stiff source and we discard all MHD fiuid back reactions onto the field itself . Prior to 
decoupling, the conductivity of the primordial plasma is very large p^ , |6^ and for practical purposes can be assumed 
infinite. In the comoving frame, this implies the "frozen- in" condition E = — v x B, where v is the plasma peculiar 
velocity and E is the electric field induced by plasma motions. Infinite conductivity leads to a vanishing electric field 
in linear perturbation theory (v <^ I) and allows the time evolution of the magnetic field to decouple from its spatial 
structure on sufficiently large scales. As the Universe expands, magnetic field lines are simply conformally diluted 
due to fiux conservation: B(?7, x) = B(x)/a^. On small scales, however, a primordial magnetic field is damped due to 
photon and neutrino viscosities |^.^ . As in Ref . [ p8| , we parametrize this damping by introducing a hard ultraviolet 
cutoff wavenumber kjj in the magnetic power spectrum. We will compute the magnetic damping cutoff wavenumbers 



feel's for both vector and tensor perturbations in Sec. Ill 



A statistically homogeneous and isotropic magnetic field must have the two-point correlation function [ p8y65| ] 

(5,(k)B*(k')) = i2TrfP,jP{k)S{k~k'), (2.1) 

where 

-Pij = ^ij kikj (2-2) 

is a projector onto the transverse plane: 

PijPjk — Piki Pij^j ^ 0' (2-3) 
and ki — ki/k. We adopt the Fourier transform convention 

B,(k) ^ I d^x exp{ik-x.)B,{x). (2.4) 



Note the projection tensor of Eq. ( pT^ ) is valid only for the case of a flat Universe where perturbations can be 
decomposed into plane waves; for non-zero spatial curvatures, the analog to a plane-wave basis must be employed 
(see, e.g., [^). A specific magnetogenesis model consists of specifying the function P{k), which we take to be a power 
law 

P{k) = Ak'\ (2.5) 

Our primary interest is to constrain the primordial comoving magnetic field strength on a certain comoving length 
scale. We therefore convolve the field with a 3D-Gaussian filter transform of comoving radius A, -Bi(k) i3i(k) * fk, 
where fk = exp(— A^A:^/2), and normalize as 

(B.(x)i3,(x))U=i?i (2.6) 

Thus B\ is the magnetic comoving mean-field amplitude obtained by smoothing over a Gaussian sphere of comoving 
radius A. The corresponding mean-square value is then given by the Fourier transform of the product of the power 
spectrum P{k) and the square of the filter transform fk, 



Bl = 7^ I d'kP{k)\M^ ^ 7:=^T;^r ( ^ ) , (2.7) 



2A 1 

(27r)3 7 " — ^■"^'•"^i - (27r)2 A"+3^ 

with the factor 2 coming from the trace of the projection tensor of Eq. ( |2.2| ). We require the spectral index n > —3 to 
prevent infrared divergence of th e inte gral over the spectrum of long wavelengths k ^ 0. Solving for the normalization 
constant A and using Eqs. ( |2.1[ ) and ( |2.5[ ), we arrive at the two-point correlation function for a primordial stochastic 
magnetic field 

(i?.(k)s;(k')) - ^ 2 r{4^f ''k^^^^~^'^' ^ ^ ^^-^^ 
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where k\ = 2tt/\. The spectrum vanishes for all scales smaller than the damping scale k > ko- The condition n > —3 
guarantees that superhorizon coherent fields are not overproduced; the limit n — s- —3 approaches a scale-invariant 
spectrum. The case n = corresponds to a white noise spectrum where we have equal power at all wavelengths. For 
a causally-generated stochastic magnetic field, we require n > 2 ||5^ , |65| , |67|] . 

The induced electromagnetic stress-energy tensor is given by the convolution of the magnetic field |pq] 



(k) 



1 



(27r)3 47r 
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Bi(p)Bj(k - p) - -5„B,(p)B,(k - p) 



(2.9) 



SB) 



n 



iS) 



n 



In the 



It can be geometrically decomposed into scalar, vector, and tensor perturbation modes, r^^ — ^^^^ 
according to their three-space coordinate transformation properties on the constant-time hypersurface 
linear approximation, all types of cosmological perturbations are decoupled from each other dynamically; thus we can 

that sources 



consider each type of perturbation independently. From the tensor n|^^ we can construct a vector n|^^ 

t(t) 



the vorticity perturbations, whereas the tensor 11-^- ' sources the gravitational wave perturbations. To obtain CMB 

derive two-point corre 
isotropic spectra |n^^^''"^^ (fc)p. This is the subject to which we now turn. 



power spectra, we need to derive two-point correlation functions for 11^ ' and 11) ' and extract their corresponding 



A. Vector Projection and Correlation Function 



We begin by illustrating how to project from a generic spatial metric perturbation Scjij its vector piece SgiJ'^ ■ A 



vector spatial metric perturbation must have the form [|69| 



where is a divergenceless three vector. A possible construction for is given by 

— k^iSg^^ii kikjYikrfidgfnn- 
The projection then follows from substituting Eq. (2.11) into Eq. ( 2.10| ): 



{km^Qmi kj^kjyikrfi5grnn)kj -\- {km^Qmj kjk^kjiSg^jui^ki — {Pinkj Pjnki)ktn^Qr) 



Using Eq. (2.12), the vector part of the electromagnetic stress-energy tensor is given by 



from which we can construct a vector II^^'' via contracting with the unit vector k 



J' 



n(^) - n^^^^- - p k t(^) 



(2.10) 
(2.11) 
(2.12) 
(2.13) 

(2.14) 



The physical meaning of 11^^^ is clear upon examining the Lorentz force vector. In the infinite conductivity limit, the 
Lorentz force vector in real space is given by |^ 

1 



L(x) 



47r 



{B(x) X [V X B(x)]} = — \ [B(x) . V] B(x) - -VB^(x) 



(2.15) 



Fourier transforming Eq. ( 2.15| ), extracting the corresponding vortical component which satisfies the divergence- 
less condition h'l^^ki — 0, and comparing with Eq. (2.14) shows that 



(2.16) 



The vector nl^-* will appear in the evolution equations for vector perturbations in Sec. IV A . 

The stochastic and transverse nature of 11^^-* lead us to define the two-point correlation function 



nr'(k)n^.^)*(k'))^p.,|n(^)(fc)p5(k~k'). 



(2.17) 
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The vector isotropic spectrum |n(^^(fc)p can be obtained using E g. (p.l4 ) for evaluating the two-poi nt co rre- 

lation function of the electromagnetic stress-energy tensor of Eq. ( |2.9| ), and comparing the result with Eq. (2.17). A 
lengthy calculation in the Appendix gives 



|n(^)(fc)|= 



1 



87r(2n -I- 3) 



i,2n+3 



k < ko- 



(2.18) 



The first term dominates when n > —3/2, whereas the second term dominates when —3 < n < —3/2. For the case 
n > —3/2, the vector isotropic spectrum becomes approximately white noise (independent of k) and depends on the 
ultraviolet cutoff wavenumber fc^. This is because the induced electromagnetic stress-energy tensor of Eq. (2_^) is 
quadratic in the stochastic magnetic field and the convolution of the magnetic field couples the large and small scale 
modes. Each mode of the vector isotropic spectrum is then affected by all scales of the magnetic power spectrum 
of Eq. (2.8) and for the case of n > —3/2, the cutoff scale perturbations completely dominate the large scale modes 
(see also Sec. V of Rcf. JM). Note that the term within the square bracket is the normalization A of the magnetic 
power spectrum in Eq. ( |2.5D . To simplify the calculation, we will only consider the corresponding dominant term for 
a given spectral index n, although including the contributions from both terms is a straightforward extension of the 
calculation presented here. In the neighborhood of n = —3/2, both terms must be included to handle correctly the 
removable singularity. 



B. Tensor Projection and Correlation Function 



Gravitational radiation is produced by the transverse and traceless piece of the electromagnetic stress-energy tensor, 
given by (see, e.g., Q) 



= fp. p. _ ip. p 



(2.19) 



It follows from the transverse and traceless properties of the tensor II,-^- ' that its two-point correlation function can 
be written as |58l 



(n(p(k)n[P*(k')) ^ M.,.„.|n(^)(fc)p5(k - k'). 



(2.20) 



The tensor structure Mijim is 

M ijlrn 



— PilPjm ~h PifYiPjl — PijPlm 

— ^il^jm + ^im^jl ^ ^ij^lm H" kikjkikm 

H~ ^ijklkjyi -f- dlmkikj — Siikjk^n ~^ ^jm^i^l ~ ^im^j^l ^jlkikj^ 



(2.21) 



and satisfies Mijij — 4 and J^uim — MijU = 0. The tensor isotropic spectrum |n(-'")(fc)p can be obtained us- 
ing Eq. ( 2.1S| ) for n|J' , evaluating the two-point correlation function of the electromagnetic stress-energy tensor of 
Eq. (2J), and comparing the result with Eq. ( 2.2C ). A similar calculation as in the case of the vector isotropic 
spectrum gives 



|n(^)(fc)|^ 



1 



167r(2n + 3) 



2T (2i±2) kl+^ 



,2n+3 



„2n+3 



k < kn- 



(2.22) 



The tensor isotropic spectrum differs from its vector counterpart only by a factor of two, due to the ratio of traces of 
their corresponding tensor structures Pij and Mijim- Again, the first term dominates when n > —3/2, whereas the 
second term dominates when —3 < n < —3/2. For the case of n > —3/2, the cutoff scale perturbations completely 
dominate the large scale modes and hence the tensor isotropic spectrum depends on the ultraviolet cutoff wavenumber 
ko'y the resulting tensor CMB temperature power spectrum then possesses the well-known behavior of a white noise 
source, I'^Ci oc P [Q. We will demonstrate that this is also true for tensor CMB polarization and temperature- 
polarization cross correlation. As above, the term within the square bracket is the normalization A of the magnetic 
power spectrum in Eq. ( |2.5| ). As with the vector case, we will only consider the dominant term of Eq. (2.22) for a 
given spectral index n. 
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III. MAGNETIC DAMPING SCALES 



The evolution and damping of primordial magnetic fields are studied in Refs. pl| j 42|. T hese authors consider cases 
for which either the magnetic field is linearized about a constant background field ^Ij or a magnetic field with a 
tangled component of unrestricted amplitude is superposed perpendicularly on a homogeneous field We are 

interested in a stochastic magnetic field with a power-law power spectrum. We will first briefly recapitulate the 
findings of Refs. |4iy42| ]. Based on these results, we then proceed to compute the magnetic damping scales separately 
for vector and tensor perturbations for a power-law magnetic field. 

Primordial magnetic fields are damped by radiative viscosity, which arises from the finite mean free paths of 
neutrinos and photons. Damping of MHD modes by neutrino viscosity is the most efficient around neutrino decou- 
pling (T ~ IMeV). At that time, the neutrino physical mean free path (Zi/dcc ~ 10^^ cm) and the Hubble length 
(iJ~dcc ~ 5 10^° cm) are comparable, hence the dissipation of magnetic energy can only occur on relatively small 
scales. Photon viscosity, on the other hand, damps MHD modes from after e+e~ annihilation (T ~ 20keV) until 
recombination (T ~ 0.25 eV); thus it is capable of dissipating magnetic energy on larger scales. There are three types 
of propagating MHD modes: fast and slow magnetosonic waves, with the fluid velocity making an arbitrary angle 
with the background magnetic field; and Alfven waves, with the fluid velocity oriented perpendicular to the wave 
vector k and the background magnetic field. Alfven waves induce neither density nor temperature perturbations. Fast 
magnetosonic waves are similar in nature to sound waves. Like the acoustically oscillating density fiuctuations, they 
are Silk-damped by radiation diffusion on scales below the radiation diffusion length. Meanwhile, slow magnetosonic 
and Alfven waves possess similar behaviors. During the radiation diffusion regime {k~y^y^ > these waves either 

oscillate negligibly or become overdamped, hence the dissipation of magnetic energy becomes inefficient. It is only 
upon entering the free-streaming regime (fcphys ^ ^1^,1) before recombination that these waves will suffer additional 
damping. The resulting maximum damping scale for these waves is dependent on the background magnetic field 
strength and is on the order of the Alfven velocity times the comoving Silk scale. Since Alfven modes describe in- 
compressible motions, we can obtain the magnetic damping scales for vector and tensor perturbations via computing 
the damping scales of such modes for a power-law magnetic power spectrum. 



A. Vector Perturbations 



Since vector perturbations induce CMB anisotropics via vorticity at recombination, we need to evaluate the Alfven 
wave damping scales at recombination for a power-law magnetic field. Around recombination, all Alfven modes are 
overdamped. The modes that suffer the most damping while overdamped are those in the free-streaming regime 
For a non- linear Alfven mode propagating in a uniform background field B, its free-streaming damping scale 
at recombination is given by Eq. (8.11) of Ref. 

where is the comoving Silk scale at recombination and we have assumed Tdcc = 0.25 eV, ft^h'^ = 0.0125, and a 
matter-dominated Universe at recombination. The Alfven velocity Va arises from the uniform background field B. 
For a linearized Alfven mode, we have to replace Va by Vacos0, where 9 is the angle between the wave vector and 
the zero-order background field (cf. Eq. (108) of Ref. pd|]). 

For a stochastic magnetic field with a power-law power spectrum, the effective homogeneous magnetic field re- 
sponsible for the Alfven velocity can be obtained via smoothing the stochastic field. As in Ref. |5^, we assume the 
field smoothed over the damping scale Xd acts as the effective homogeneous field i?cff. For each spectral index n, 
Bl cx A-("+3) [cf. Eq. {^)], and B^g is related to Bx through (see also Eq. (26) of Ref. ||) 

Tl+3 

B^s^bJ^) ' . (3.2) 



Smoothing the stochastic magnetic field on scales larger than the damping scale will result in a smaller effective 
homogeneous field, hence a smaller effective Alfven velocity and a larger momentum cutoff wavcnumber kjj. Since 
for n > — 3/2 the vorticity source becomes approximately white noise (independent of fc) and is fcu-dependent [cf. 



Eq. (2.18)], a larger momentum cut off wavenumber ko will give rise to larger CMB an isot ropics in t his regime as we 



will see. The estimation in Eq. (^) is therefore a conservative one. Substituting Eq. (|3.2| ) into (|3.l|) yields 
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Note that for a given spectral index n, Bf/k'^' 



in Eq. 



cx A, where A is the normahzation of the magnetic power spectrum 



B. Tensor Perturbations 



Since the sourcing of gravitational radiation after the Universe becomes matter-dominated is negligible (see Sec. IV B 
and also Ref. |5^), the relevant tensor damping scales are the Alfven wave damping scales at equality. As in recom- 
bination, all Alfven modes are also overdamped around equality, hence the modes that are undergoing free streaming 



suffer the most damping. The situation is clearly depicted in Fig. 1 of Ref. 
damping scale at equality is (cf. Eq. (106) of Ref. Q and Eq. (8.10) of Ref. 



1,-1 _ :1R 



The Alfven wave free-streaming 



(3.4) 



where L^'*(Teq) is the photon comoving diffusion length at equality. 



L^'S(Teq) « 19.5 



cq 



-5/4 



-1/2 



0.0125 



assuming T^q = 5.5 eV(rio^^), — li and Qi,h^ 
as in the case of vector perturbations gives 



0.25 eV 

0.0125. Substituting Eq 



-1/2 



Mpc 



(3.5) 



into (3.4), a similar manipulation 



.3 X 10-^)" 



10-9 G 



kx 



IMpc" 



1 + 3 
1 + 5 



Mpc 



(3.6) 



Note again that for a given spectral index n, B^/k^ 
spectrum in Eq. (2.5). 



2/I.n+3 ^ A 



where A is the normalization of the magnetic power 



Several comments are in order. First, since the tensor source contributes earlier (at equality) than the vector source 
(at recombination), tensor perturbations are damped at smaller scales as illustrated by Eqs. (3.3) and (3.6). For 
n > —3/2, CMB power spectra are dependent on the momentum cutoff wavenumber fc^i and scale with it as fc^"'"'^, 
we therefore expect tensor perturbations to generate larger anisotropics than the vector perturbations in this regime, 
at least for I < 500 that we are considering. Second, in Ref. ||5^, the magnetic damping cutoff wavenumber for tensor 
perturbations is found to be 4.5 Mpc^^. This value is derived based on the assumption that the Alfven modes are 
undergoing damped oscillatory motions. Our analysis however, shows that for the magnetic field strengths considered 
here, the Alfven modes should be in the overdamped free-streaming regime around equality, as also illustrated in 
Fig. 1 of Ref. I^lj] . Finally, as poi nted out in Ref. |4^ , the Alfven damping scale at equality could in principle be 
larger than that given by Eq. (3.4) since additional damping could arise due to a possible breakdown of the WKB 
approximation in the regime where the Alfven mode is undergoing overdamped fr ee s treaming. In the absence of 
an accurate quantitative treatment for Alfven damping scales in this regime, Eq. ( [3.4| ) is our best-educated guess. 
Nevertheless, we caution the readers that with a possible larger damping scale fc^^ than that given by Eq. (3.4), the 
induced tensor anisotropics for n > —3/2 will be reduced accordingly. 



IV. METRIC PERTURBATIONS AND THEIR EVOLUTION 



A primordial stochastic magnetic field generates CMB anisotropics via its gravitational effects on the metric tensor. 
The full metric tensor can be decomposed into its background and perturbation pieces, 5^,^ — gj^J + Sg^i,; for a flat 



Universe with the usual conformal FRW metric, 



(0) 



where ry^^ = diag(— 1, 1, 1, 1) is the Minkowski metric 

tensor. The vector (Sec. IV A ) and tensor (Sec. IV B| ) perturbations are calculated separately; scalar perturbations 



will generally result in smaller CMB anisotropics compared to vector and tensor contributions, as argued in Sec. VIII 
and so will not be considered here. We review the various metric tensor contributions and give the corresponding 
evolution equations due to a primordial stochastic magnetic field. We then obtain solutions to these equations, which 
can be expressed as functions of the isotropic spectra derived in Sec. |l|. 



7 



A. Vector Perturbations 



Vector pert urbat ions to the geometry are described by two divergenceless three-vectors d and ^ with the general 
form [see Eq. ( 2.10| ) also] 



Sg'ol ' = Sg^p = a^^^k, + S^fh). (4.1) 

Vector perturbations represent vorticity; the divergenceless condition for vectors and guarantees the absence of 
density perturbations. Vector perturbations exhibit gauge freedom, which arises because the mapping of coordinates 
between the perturbed physical manifold and the background is not unique. From vectors Q and ^i, we can construct 
a gauge-invariant vector potential Vi = Ci~^d/f^ that geometrically describes the vector perturbations of the extrinsic 
curvature | |70| , |7l[| . We now exploit the gauge freedom by explicitly choosing to be a constant vector in time; it 

follows that Sg^'^ = —a?Vi. Vector perturbations of the stress-energy tensor can be parametrized by a divergenceless 
three-vector v^^^ that perturbs the four-velocity = (1, 0, 0, 0) of a stationary fluid element in the comoving frame 

5u^ = {Qy'^la). (4.2) 

We can now construct a gauge-invariant, divergenceless three- vector termed the "vorticity," 

n,^vf^-V,. (4.3) 

Two Einstein equations govern vector perturbation evolution. The first describes the vector potential evolution 
under the influence of a primordial stochastic magnetic field: 

y.(,,k) + 2^F.(,,k)=- ^^-^^f(^) , (4.4) 



where n^^-'(k) is given by Eq. (2.14) and we neglect the vector anisotropic stress of the plasma, which is in general 



negligible. The magnetic field source terms (k) and 11^^^ (k) are expressed in terms of present comoving magnetic 
field amplitudes. Since both of these terms depend on the magnetic field quadratically, the explicit time dependence 
of the magnetic stress is given by n(?7, k) = n(k)/a'*. In the absence of the magnetic source term, the homogeneous 
solution of this equation behaves like Vi cx 1/a^. The complete solution including the magnetic source is simply 

During the radiation-dominated epoch we have a oc 77; a magnetic field therefore causes vector perturbations to decay 
less rapidly (1/a instead of 1/a^) with the Universe's expansion. The second vector Einstein equation is a constraint 
that relates the vector potential to the vorticity: 

- fcV,(77,k) = 167rGa\p + p)n,{7j,k). (4.6) 

Vector conservation equations can be obtained via covariant conservation of the stress-energy tensor. Since vector 
perturbations cannot generate density perturbations, we have 

d^=db = 0. (4.7) 

Before decoupling, photons are coupled to baryons via Thomson scattering. The magnetic field affects the photon- 
baryon fluid dynamics via the baryons; wc therefore introduce the Lorentz force term into the baryon Euler equation. 
The Euler equations for photons and baryons are respectively |^,^ 

fi7»+T(z;^r-«ir)=o, (4.8) 

+ -^u - Uv^' - v'^') = f'^y (4.9) 

In the above, fi^^b — v^^'' —V represent vorticities of photons and baryons; r = n^^aTa is the differential optical depth 
where rie is the free electron density and ctt is the Thomson cross section; R = (pt +Pb)/{p-f + P-y) — 3pfc/4p^ is the 
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mom entu m density ratio between baryons and photons; and is the vortical piece of the Lorentz force given by 



Eq. (2.16). Again, we neglect the small effects due to the vector anisotropic stress of the plasma. This set of vector 



conservation equations are similar to the one that describes Alfven waves in Ref. Equations (14), (f^), ([f.q), 

and ( |4.9| ) are not independent. Using the definitions of i?, L^^^\ and the fact that [p^ +P7) « and solving the 

Euler equations in the tight-coupling approximation v^^^^ ~ ^ obtain the following approximate solution for 
the vorticity: 

Note that the same result can be obtained using Eqs. (j^^) and ([4.6[ ). The factor 1 + R represents reduction in 
the vorticity due to the Compton drag of baryons. At decoupling, the momentum density ratio between baryons 
and photons has an approximate value of i?doc — 3/3bo/4p7oZdcc — 0.35, where we have assumed z^icc = 1100 and 
flbh^ — 0.0125. The vorticity solution of Eq. ( 4.10| ) is valid for perturbation wavelengths larger than the comoving 



Silk scale Ls, where photon viscosity can be neglected compared to the Lorentz force. For k > ks, where kg = 27r /Ls, 
the Euler equation that includes the viscous effect of photons is 

-p^ + pA + I p,- + ^) ^ ^ y , (4.11) 



where x = (4/15)p^-L^a is the photon shear viscosity cofhcient and = is the photon comoving mean-free 
path. In this regime, the vorticity can be obtained using the terminal-velocity approximation. Equating the photon 
viscosity term to the Lorentz force, we obtain Q 

k) c I ^ ' -, k > ks. (4.12) 

The next step is to introduce two-point correlation functions for the vector potential and the vorticity. Defining 
their two-point correlation functions as in Eq. ( ^.17 ) for the vector n|^\ and taking ensemble averages of Eqs. (|4.5|), 



(4.10), and (4.12) the rms isotropic spectra for the vector potential and the vorticity are simply 



fcn<'^'(fc)?7 



k < ks 



niv,k)^i ^^-^w.o^P.o; _ _ (4.14) 



k > ks- 



{kLy/5)(p~,o+p~,o) 

Vector perturbations induce CMB temperature anisotropics via a Doppler and an integrated Sachs- Wolfe effect : 

rno 

e(^)(7?o,k,n) = -v(^) -ni;;"^^ + / drjV-h, (4.15) 

where ?7dcc represents the conformal time at decoupling. The decaying nature of the vector potential V implies 
that most of its contributions toward the integrated Sachs- Wolfe term are around T^dcc- Neglecting a possible dipole 
contribution due to v^^^ today, we obtain Q 

e(^) (770, k, n) ~ v(^) (?7dcc, k) • n - V(r7dec, k) • n = 0(77dec, k) • n. (4.16) 
Vector CMB temperature anisotropics are due to the vorticity at decoupling. 

B. Tensor Perturbations 

Tensor perturbations to the geometry are described by 

Sg'^P =2a^h,j, (4.17) 
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where hij is a symmetric, transverse {hijkj = 0), and traceless {ha = 0) three-tensor. Unhke vector perturbations, 
tensor perturbations have no gauge freedom. 

The tensor Einstein equation that describes the evolution of gravitational waves sourced by a stochastic magnetic 
field is 



/i,,(7?,k) + 2-/i,,(ry,k) + k^h,,{rj,k) = 8nGn^\k)/a^, 
a ■' 



(4.18) 



(T) n k 

where 11^ (k) is given by Eq. (2.19) and as in the case of the vector perturbations, we neglect the tensor anisotropic 



stress of the plasma, which is in general negligible. Gravitational waves induce CMB temperature anisotropics by 
causing photons to propagate along perturbed geodesies : 



i-Vo 

e'^)(r/o,k, n) ~ / dr]hij{r],'k)ninj 



(4.19) 



Our task is, therefore, to obtain the solution for hij. To calculate tensor CMB power spectra, we need to define 



two-point correlation fmictions for hij and hij as in Eq. ( 2.20| ) for the tensor n|J\ with rms isotropic spectra h 

and h respectively Solutions to the homogeneous equation with n(-^^(fc) — are easily obtained. During the 
radiation-dominated epoch, a cx 77 and h cx jo{krj) or yo{kr]), while during the matter-dominated epoch, a oc 77^ 
and h cx ji{krj)/krj or yi{k'q)/k'q^ where ji and yi are the usual spherical Bessel functions. Assuming the primordial 
stochastic magnetic field is generated at yyin, a Green function technique yields the following inhomogeneous solution 
for the radiation-dominated epoch: 



27rGn(^)(fc); 



2 2 
cq^cq 



(3 - 2V2)k 



V 



, sin[fc(?7 
dr] — 



r]' 



(4.20) 



where rjcq denotes the conformal time at matter-radiation equality. The magnetic source term on the right hand 
side of Eq. (4.18) decays more rapidly with rj in the matter-dominated epoch than in the radiation-dominated epoch. 
An approximate solution, therefore, can be obtained by matching the radiation-dominated inhomogeneous solution 
of Eq. ( 4.2OD to the matter-dominated homogeneous solutions at equality. Retaining the dominant contribution, we 
obtain [dST 



hirj,k) ~ 47rG?7o'^cqln ( ^) fcn(^)(fc)^4^ 



krj 



(4.21) 



V. TEMPERATURE POWER SPECTRA 

We employ the total angular momentum representation introduced by Hu and White ||5^ to compute the CMB 
power spectra induced by a primordial stochastic magnetic field. By combining intrinsic angular structure with that of 
the plane- wave spatial dependence, this representation renders a transparent description of CMB anisotropy formation 
as each moment corresponds directly to an observable angular sky pattern via its integral solution of the Boltzmann 
equations. The CMB temperature power spectrum today is given by Eq. (56) of Ref. pfl]: 

ee™^lf,,,,e!^e!^^ ,5.1) 



where X stands for V or T, and Qi's are the temperature fiuctuation AT/T moments. Note that Eq. ( |3.lD is larger 
than the corresponding expression in Ref. by a factor of two as we have already taken into account the fact that 
both vector and tensor perturbations stimulate two modes individually, corresponding to m = ±1,±2 respectively in 
the notation of Ref. M]. Our strategy is to evaluate the Boltzmann temperature integral solutions to obtain the Oj's 



due to the vector and tensor perturbations. We then substitute them into Eq. (5.1) to yield the corresponding CMB 
temperature fluctuations spectra. Though the tensor results are already given in Ref. p8[ , the vector results derived 
here are new. 
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A. Vector Temperature Power Spectra 



The Boltzmann temperature integral solution for veetor perturbations is given by Eqs. (61) and (74) of Ref. 



21 + 1 



r 



fp(^)jf^)[fc(,?o-r])]}, 



(5.2) 



where 



p(V) 



^/ik 



(5.3) 



is the vector polarization source that is generated when tight coupling breaks down on small scales, where the photon 
diffusion l engt h and the perturbation wavelength become comparab le. T he approximation in Eq. (5.3) is obtained 
using Eq. (4.3) and noting that Vl dominates V at decoupling [cf. Eqs. (4.13) and (4.14)] for k > 0.006 Mpc~^, resulting 
in V co ntrib uting negligibly compared to 17 upon integrating over /c's to obtain the vector temperature power spectrum 
in Eq. (5.1). Unlike scalar perturbations, vector perturbations cannot produce compressional modes due to the lack 
of pressure support. In the usual case of vector perturbations in cosmological fluids without a magnetic field, tight- 
coupling expansion of photon and baryon Euler equations give « V, resulting in the vector polarization source 
being dependent on the vector potential instead (see Eq. (94) of Ref. [£9[). A primordial stochastic magnetic field 
thus enhances vector polarization by sourcing the vorticity. The vector temperature radial functions jp^' and 
which describe how distant sources contribute, are given by Eq. (15) of Ref. ||59[|: 



l{l + l)3i{x) 



^.(2y)^^^ ^ _/3Z(Z + l) d fji{x) 



(5.4) 



2 dx 

The optical depth between rj and ryo is defined as T{r]) = J^" drj' t{i]'), thus dr/dij — — f. Integrating Eq. (5^) by 
parts using de^'^/dri = fe^'^ and jl^^\x) — A/3(j}^^'(x))' and Eqs. ( ^ ) and ( |3.3[ ), we obtain 



21 + 1 



9 f 



(5.5) 



For the usual vector perturbations without a magnetic field, the term proportional to J;*"^^' is strongly suppressed since 



V at decoupling M] and hence f2 ~ as mentioned above. Here we have a primordial stochastic magnetic field 



sourcing fl; the term proportional to jf^^'' is then suppressed relative to the term proportional to j^''^ ' due to the 
factor fc/f. Moreover, jf^^'' has less angular power compared to j^'^'^^ (see Fig. 3 of Ref. ||5^). Thus to simplify the 

calculation, we consider only the term proportional to j;^^^^ in computing the vector temperature integral solution. 
Including the small corrections due to the angular dependence of polarization coming from the term proportional to 
P'^^jp^-* and also the vector potential will yield an additional contribution of at most a few percent toward our final 
estimate of the vector temperature power spectra. The combination fe""^ is the conformal visibility function, which 
represents the probability that a photon last scattered within drj of 77 and hence is sharply peaked at the decoupling 
period. For I < 500, we can approximate the vector temperature integral solution reasonably well as 



(IV) 



21 + 1 



ji{kr]o) 
km 



(5.6) 



using Eq. (5^) and the fact that rjo ^ JTdcc- The vector CMB t emp erature anisotropics are due to the vorticity 
at decoupling, as also illustrated by E n. (]4.16| ). S ubsti tuting Eq. (5.6) into the CMB temperature power spectrum 
expression of Eq. ( |5.l| ) and using Eqs. ( |4.14[ ) a!nd ( ^.18 ), we obtain 



C 



ee{v} 
I 



(2^)2 



-1(1 + 1) 



■"AX 



{kD-qo) 



2n+3 



r2 »±i (2n 



k 

ko 



2 

2n+3' 



72 



3) (fcAr?o)2"+6 

i/2ikm). 



(1 + i?dcc)^ 



fes 



dk k 



25 

7 dec 



fes 



dk 



(5.7) 
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where we have defined the Alfven velocity as vax = i?A/[47r(p^o +^70)]^^^ — 3.8 x 10~^(i?A/10~^ G). Note that 
(271)^"^^^ v\-^^/ [r^ (^^^) ^a"^^] where A is the normaUzation of the magnetic power spectrum in Eq. (2.5). 

Depending on whether n > —3/2 or —3 < n < —3/2, we retain only the corresponding dominant term of the vector 
isotropic spectrum in Eq. ( |5.7| ). First consider the case n > —3/2, where the vorticity source becomes approximately 
white noise (independent of k) and that |n(^^(fc)p is dependent on k^- To obtain an analytic estimate of the integral 
fl^'^ dkkJi_^^ ^2(^^0)1 consider the more general integral J^^ dxxPJf{x) for some p > 0. Since Jf{x) only begins to 
contribute to the integral significantly when a; ^ Z, in this limit, we employ in the integral the Ji{x) asymptotic 
expansion for large argument |72j: Ji{x) ^ ^J2/ (ttx) cos[a; — {21 + l)7r/4]. Approximating the oscillations by a factor 
of one half then gives 

dxxPJ^(x)~ dxxPJ^(x) ~ I / ^ (5.8) 

The approximation tends to underestimate; it is good to a few percent for p > 1 and is within 30% for < p < 1. 
The integral dkkJfj^^^^.i^'Ho) corresponds to the case p = 1; the remaining integral on the rhs of Eq. (5.7) 

dk k~^ J^j^^i^ikrio) can be well-approximated by {k^'^ — fc^"^)/ (371770), which is good to within 20%. As in Eq. (^.^), 
this approximation is obtained via employing the Bessel function asymptotic expansion for large argument, which is 
justified since fcr^o > I for ks < k < ko and / < 500. Keeping only the highest-order term in Z, we obtain the vector 
CMB temperature power spectrum for n > —3/2: 

^2(^ee(v) ^ (27r) v^^l (kDm) 

I o 



r2(ii±3)(2n + 3) (fcAr/o)2"+6 




25 / f]o 



1 1 



3 V-^7dcc/ [{ksrioy {komY 



n>-3/2. (5.9) 



The dominant contribution comes from the term proportional to {ksijo — l), which arises from the non-damped vorticity 
for k < ks- The remaining term arisin g fro m the damped vorticity for ks < k < kjj gives a negligible contribution of 
< 1%. A numerical evaluation of Eq. (5.7) shows that its second integral on the rhs, which arises from the damped 
vorticity, always contribute negligibly (< 1%) compared to its first integral for I < 500 and all cases of n. Therefore, 
we will neglect the damped vorticity contribution when evaluating the remaining vector CMB temperature power 
spectra. 

For —3 < n < —3/2, the needed integral is J^^^ dfc fc^"+'* Jj^^j^^2(^'7o)- We must consider three cases depending on 
whether the exponent 2n -I- 4 is greater than, equal to, or less than zero. For —2 < n < —3/2, using Eq. ( ^.8[ ) for 
p = 2n -|- 4, we obtain 

,2^ee(V) _ (2^)'"+'" f Vdcc/VO V Vax^I^ r,, „ x2„+4 _ ,2n+41 

' 4 l^l + i?decy' r2(2±3)(2n + 3)(n + 2)(n + 3)(fcA%)2"+6 ^' 

-2<n< -3/2. (5.10) 

For n — —2, again using Eq. ( pT^ ) for p = 0, we obtain 

fcr^^^ = 2(2.)^ ( ^!^) ' ^ In f ^ V n = -2. (5.11) 

For —3 < n < —2, numerical evaluation shows that for —2.3 ^ n < —2, the integral J^^^ c?fc J;-^^j^^2(^%) can be 

well-approximated by [{ksrjo)^^^^^ — Z^"+^]/[27r(n-|- 2)77^"^'^], which underestimates by at most 30-40%. The resulting 
temperature power spectrum is then formally identical to that of the case —2 < n < —3/2: 

,2pee(v) ^ (27r)^"+i° / rjdec/vo V Vax^I^ r., „ N2n+4 _ ,2«+4] 

' 4 U + i?decy r2(^)(2" + 3)(n + 2)(n + 3)(fcAr?o)2"+6 J' 

-2.3<n<-2. (5.12) 

For —3 < 71 < —2.3, the dominant contribution to the integral J^^ dkk^'^^'^jf^^^^ihrjo) is coming from long wave- 
lengths fc — > 0, we therefore approximate by integrating over k to infinity. The resulting integral can be evaluated 
analytically by using 6.574.2 of Ref. 171, 
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,6-1 



dk Jp{ak)Jg{ak)k = 



T{b)T 



p+q-b+l 
2 



and 8.335.1 of Ref. 



r(2a;) 



-r{x)r [x + - 



Keeping only the highest-order term in /, we finally obtain 



r2(-n- 



Re (p + q + 1) > Re 6 > 0, a > 0; 

(5.13) 

(5.14) 

-3<?i<-2.3. (5.15) 



22n+7 yi + R^^J r(-2n-4)r2(li±3) (2n + 3)(n + 3)(A:Ar;o)2"+6' 

Our approximation overestimates, as expected, and the accuracy improves as n decreases since more contribution 
arises from small k and hence the result will be less sensitive to the upper limit of the integral, which we have 
approximated to be infinity. It is good to within 30% for —2.5 < n < —2.3 and a few percent for —3 < n < —2.5. 
The temperature power spectrum of each case above has the same I and fc^i (with k\,ks k^) dependence as 
the corresponding spectrum induced by a primordial homogeneous magnetic field p6| (the correspondence between 
the spectral index of Ref. |4^] and ours is n — > 2n + 3). For the sake of completeness, we note that the vector 
potential contribution arising from the jp^'' term in Eq. (^.5D will induce temperature power spectra Pcf ^^^^ oc 



for n > -3/2 and PC, 



2^ee(y) 



oc 



We now consider the case n = —3/2 a nd sh ow that this apparent singularity is removable by considering both terms 
of the vector isotropic spectrum in Eq. (2.18). In the limit n = —3/2 + e, we have 



1 



16 r2(3/4)/c3 2e 



, 2eVV 



2e 



(5.16) 



Upon expanding the expression within the square bracket to 0(e) and using the small- a; expansion to the first order, 
i.e. ln(l + x) ~ X for x = {k — k]j)/kD, we obtain 



|n(^)(A)|^ 



{2n) 



16 r2(3/4)fc3 V3 



k 

ko 



-3/2. 



(5.17) 



The same result can be obtained via direct substitution of n = —3/2 in Eqs. (A9) to (All). Using Eqs. ( ^.17 ) and 
(WM) for p = 1 and 2, a similar calculation as in Eq. (|5.9|) gives 



p^ee(v) ^ (27r) f rjdec/vo 



l + i?do 



r2(3/4)(fcA?7o)3 



10,, {ksm? - 

3 korjo 



-3/2, 



(5.18) 



thus showing that the singularity at ti = —3/2 is indeed removable. For the rest of the paper, we will not pro- 
duce explicit power spectrum expressions for the case n = —3/2. Readers who are interested can easily derive the 
corresponding results via a straightforward extension of the calculation outlined above. 



B. Tensor Temperature Power Spectra 

The Boltzmann temperature integral solution for tensor perturbations is given by Eqs. (61) and (74) of Ref. 



21 + 1 



where 



piT) ^ 

3f 



(5.19) 



(5.20) 
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is the tensor polarization source and jp"^^ is the tensor temperature radial function given by Eq. (15) of Ref. p9| 



Jl (X) 



l3 {l + 2y.ji{x) 
8(/-2)! a;2 ' 



(5.21) 



Using Eq. (4.21) and defining x = krj and xq = fc?7o, we approximate the tensor temperature integral solution as (see 
also Eq. (ISjof Ref. ||) 



er(>7o,fc) 

21 + 1 



, 80 + 2)! 



nm(fc) /'^°rf^^2(x)jKxo-g_ 

X [xq — xY 



(5.22) 



The integral above can be numerically approximated as in Eq. (19) of Ref. 



dx 



j2{x) jljxQ ~ X) _ TT 

X {xq — xY 2 
7tt 



dx 



■h/2{.x) Jl+1/2{XQ - X) 

a;'^/^ {xq — a;)^/^ 



— Vl / dx 
20 



-^5/2(2;) J/+i/2(a:o - x) 



X {xq — xY 



Ttt [31Ji+3{xq) 



50 V 2 xl 



(5.23) 



where in going from the second to the third line, we have inserted a factor of 1/3/2 for better numerical agreement 
and used 6.581.2 of Ref. f^: 



f 

Jo 



dxx^ ^{a — x) ^ Jp{x)Jq{a — x) 



(-l)"r(6 + P + m)T{b + m) 



ryb °0 
- V 

aq ^ m!r(6)r(» + m + 1) 

m— \ / \-t / 

Re(6 + p) > 0, Reg > 0. 



(6 + j9 + g + 2m) Jfc+p+,+2„i(a), 



(5.24) 



Numerical evalua tion shows that the approximation in the second line of Eq. ( 5.2g ) is good to 10% for I < 500. 
Substituting Eq. (|5.23|) into Eq. (|5.22|) yields 



2^ + 1 



(l^2)\ 



G77oZoqln — 



n(^)(fc) 



{kf]oY 



(5.25) 



Using Eqs. (|2.22[) , (jsll ), and (|5.25| ), we obtain 



ee(T) 
I 



49 
10000 

X 



(27r)^"+^-'P(/- l)(/ + l)(; + 2) 







n + 3 V fc£) 



2ri+3 



Grjlzeq In 



cq 



n 2 



(kDm) 



2n+3 



r2(ii±3)(2n + 3)77i5 (fcAr?o)2"+6 



(5.26) 



Note that (27r)^"+^"i3|/ [F^ (^^) ^^"^^] oc A^, where A is the normalization of the magnetic power spectrum in 
Eq. (U). 

For n > —3/2, the gravitational wave source is fcu-dependent, and the resulting temperature fluctuation spectrum 
possesses the well-known behavior PCi oc l^. The integral J^" dkk~'^ jf^^{kriQ) can be evaluated using Eq. ( 5.13 ); 
retaining only the highest-order term in Z, we obtain 



^2(^ee(T) 



49 
7500 



(2^) 



2TI+11 



GriQZcq In 



(koVo) 



2n+3 



r2(«±3)(2n + 3) {kxT]oY"+e 



n > -3/2. 



(5.27) 



For —3 < n < —3/2, a similar calculation gives 



^2^ee(T) ^ 22"-5 — (2^)2"+i2 
' 625^ 



Gr]QZcq In 



F(l - 2n) 



Bin 



r2(l - n)F2 (li±3) (2n + 3)(n + 3) \kxTjo 



2n+6 



-3 < ?i< -3/2. 

Equivalent tensor perturbation results are given in Eqs. (20) to (22) of Ref. [sl 



(5.28) 
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VI. POLARIZATION POWER SPECTRA 



Polarization of the CMB comes in two flavors: E-type and B-type with electric (—1)' and magnetic (—1)'+^ parities 
respectively [^,^. Physically, they represent polarization patterns rotated by 7r/4 due to the interchanging of Q 
and U Stokes parameters. Vector and tensor perturbations induce both types of polarizations. Scalar perturbations, 
however, cannot generate B-type polarization due to azimuthal symmetry. A detection of the B-type polarization 
from future high sensitivity CMB polarization measurements therefore would provide compelling evidence for vector 



and/or tensor contributions. Similar to the CMB temperature power spectrum of Eq. (5.1), the E-type and B-type 
polarization power spectra are respectively 



' TT y 21 + 1 21+1 ' ^ ' 

where X stands for V or T. The correspondence between notations of Ref. |5|] and ours for polarization moments 
are Ef^'^'^ E^^^ and B^^^'' ±Bi^\ and similarly for the tensor perturbations. 



A. Vector Polarization Power Spectra 



1. E-type Polarization 
The E-type polarization integral solution for vector perturbations is pE 



= -Vel" d,fe-^pi-hP[ki,o - V)], 



where 



(x) = -v/a-i)(z + 2) 



(6.3) 



(6.4) 



is the vector E-type polarization radial function given by Eq. (17) of Ref. Q. Using Eq. (U) for F(^) and the 
spherical Bessel function recurrence relation 



- j'lix) =ji+i{x), 

we approximate the vector E-type polarization integral solution as in Eq. ( |3.6| 

El''\r^o,k) ,/(/-l)(/ + 2) 



21 + 1 



18 



-kLj dcc^('7doci k) 



jiikyq) _ ji+i{km) 



km 



(6.5) 



(6.6) 



where L-ydcc — T'doc — 3-39 Mpc is the photon comoving mean-free path at decoupling, assuming Tjcc — 0.2 5 eV and 
flbh^ — 0.0125. For t he tig ht-co upling approximation to be valid, we require kLj^cc < 1- Substituting Eq. (|6.6| ) into 
(3.1) and using Eqs. (4.14) and (2.18), we obtain 



-Rdc 



■"AX 



(koVo) 



2n+3 



7 dec 



ks 



dkk^ 



n + 3 \k 



r2(li±3)(2n-f 3) (fcA??o)2"+6 

. Jl+l/2{kTJo) Ji+3/2{kvo) 
^ ' [km? km 



(6.7) 



As in the computation of the vector temperature power spectra in Sec. V A , we have neglected the damped 
vorticity term, which again contributes negligibly (< 3%) for I < 500 and all cases of n. No te t hat 
{2tt)'^^^^°v\-^/ [r^ (^^) fc^"^^] oc A'^, where A is the normalization of the magnetic power spectrum in Eq. (2.5). 
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Again, depending on whether n > —3/ 2 or —3 < n < —3/2, we retain only the corresponding dominant term of 
the vector isotropic spectrum in Eq. ( |6.7|) . A further simphfication occurs by noting that although the cross term 
proportional to o//+i/2(^?yo)^;+3/2(fc'?o) is difficult to evaluate analytically, a numerical evaluation shows that its value 
is approximately minus twice that of the term proportional to ^;+i/2(^'?o) cases of n. First consider the case 

n > —3/2, where the vorticity source becomes approximately white noise and is fc^i-dependent. The relevant integrals 
are 



{l + lf dxxjf+y^{x) 



and 



(6.8) 



where we have defined x = kijo and xs = ksrjo- Using Eq. ( p.SD for p = 1 and 3 respectively for these two integrals, 
we obtain 



2r<EE{V) _ (27r) 



2n+10 



dec 



18 V 1 + ^dc 
n > -3/2. 



2n+3 



r2 (li±3) (2n + 3) (fcAr?o)2"+6 



i\ksm-i) 



(6.9) 



Comparing to the numerical evaluation of Eq. ( p.7[ ) shows that our approximation is good to a few percent. 
For —3 < n < —3/2, we need to evaluate 



dx x^"^^"^ Jfj^^ I2{x) 



and 



(6.10) 



Since the exponent within the first integral 2n + 4 changes sign whereas the exponent within the second integral 2n + 6 
remains positive throughout —3 < n < —3/2, as in the vector temperature power spectra calculation, we consider 
cases depending on whet her 2 n + 4 is greater than, equal to, o r le ss than zero in this regime. For —2 < n < —3/2, 
the two integrals of Eq. ( 6.10| ) can be approximated using Eq. ( |5.8| ) for p = 2n + 4 and 2n + 6 respectively, hence 



l2cEE(V) ^ (27r)^"+l" ( 77dec/r/0 



36 



l + i?dc 



2 / r \ 2 

dec 



2n+6 



(fcs7yo)2"+6 - ;2"+6 (fc5??o)'"+^ - 



-2 < n < -3/2. 



(6.11) 



Here our approx imat ion is good to within 10%. For n 
integrals of Eq. ( 3.1C| ) , we obtain 



2, using Eq. (5.8) for p ~ and 2 respectively for the two 



27r)5 / ridec/vo 



9 



l + i?do 



-^7 dec 



""AX' 



[kxrioY 



{ksmf 



(6.12) 



Comparing to the numerical evaluation of Eq. ( |6.7| ), the approximation here is good to within 10 % in general. For 
—3 < n < —2, as discussed in Sec. V A , a numerical evaluation shows that th e firs t integral of Eq. ( |6.10| ) can be well- 
approximated by P[{ksr]Q)^"'+'^-P"^/i2n+'i)T: for -2.3 < n < -2 [cf. Eq. ( ^.12D ] whereas for -3 < n < -2.3, it can 
be approximated using Eq. (5.13) since the dominant contribution to the integral arises from long wavelengths fc — > 
[cf. Eq. ( 5.15D ]. The approximation of the second integral of Eq. ( |6.10| ) usi ng E q. ( |5.8D tends to underestimate. This, 
however, can be compensated via approximating the first integral of Eq. (3.1C) by / [(fcs'7o)^"^'' ^ P"^^]/{2n + 4)7r 



throughout the regime —3 < n < —2. The resulting vector E-type polarization power spectrum is then formally 
identical to that of the case —2 < n < —3/2, with accuracy good to within 15% in general. Hence 



2qEE{v) _ (2vr 



,2n+10 



??dcc/??0 

36 Vl + ^dcc 



L 



7 dec 



-I 



I (ksVo) 



p2 ^7l+3^ 

2n+4 _ pn^ 



Vo ) (^j (2n + 3)(n + 3)(fcA?7o) 



n - 



-3 < n < -2. 



(6.13) 
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2. B-type Polarization 
The B-type polarization integral solution for vector perturbations is [ p9[ 

where 



^W(^) = _!)(/ + 2) 



(6.14) 



(6.15) 



is the vector B-type polarization radial function. Using the same approximation in Eq. ( |6.14 ) as in Eq. (5.6), we 
obtain 



21 + 1 "V 18 -"ru..-v,>..-.vv , 

which upon substituting into Eq. (6.2) and using Eqs. ( pl4| ) and ( pl8| ), yields 



(6.16) 



BB{V) 
I 



36 



X ^7 dec 



-a-i)a + 2) 



1 + -Rdct 



i2n+3 



3 V fc 



r2(«±3) (2n + 3) (fcAr/o)2"+6 

2ri+3' 



'^f+l/2(fe^0) 

(fc77o)2 



(6.17) 



where again we have neglected the contribution coming from the damped vorticity term, which is negligible (< 4%) 



for / < 500 and all cases of n. No te t hat {2tt) 



2n+10,,4 



/[r^ {^^) oc A^, where A is the normahzation of the 



magnetic power spectrum in Eq. ( p.5|) . Except for the order of the Bessel function, Eq. ( 6.17 ) is i den tical to the term 



proportional to ^1^3/2(^^0) in the vector E-type polarization power spectrum expression of Eq. (S.7). For n > —3/2 



using Eq. ( |5.8D for p = 3, we obtain [cf. Eq. (|6.9|)] 



2r<BB{v) _ (27r 



i2n+10 



54 



1 + i?do. 



For —3 < n < —3/2, the exponent within the integral 2n -f 6 remains positive throughout; thus using Eq. (5 
p = 2n 4- 6, we obtain [cf. Eq. (3.11)] 



for 



2^SB 



I'C, 



(V) _ (2^)^"+i° / 



36 Vl + ^dcc 
-3 < n < -3/2. 



2 / T \ 2 

-^^7 dec 



'70 



4 74 



p2 (n±3) (2n -f 3)(n 3)2(fcA??o) 



2n-|-6 



[(fcs% 



,2n+6 72n+61 



(6.19) 



Our accuracy here is good to the quality of the analytic approximation in Eq. (5_^) and is always within 20%. 



B. Tensor Polarization Power Spectra 

1. E-type Polarization 
The E-type polarization integral solution for tensor perturbations is [59|] 



21- 



^ Jo 



where p(^) is given by Eq. ( 5.20 ), and 



-ji{x)+j'/{x) + 2 



(6.20) 



(6.21) 
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is the tensor E-type polarization radial function. Using Eqs. (4.21) and (6.5) and the spherical Bessel function 
recurrence relation [ [72| 



l + l 



jiix)+j'i{x) =ji-i{x), 



(6.22) 



and defining x = krj and xq = fcryo, we approximate the tensor E-type polarization integral solution as 

27r 
76 



EP{r,o,k) 2tt 



21 + 1 



GryoZcqln ( — 



Jo 



X 



eq 

-2 



n(^)(/fc) 
(^ + i)(^ + 2) 

(2:0 - x)^ 



ji{xo - x) 



Xq — X 



-ji+iixo - x) 



A similar manipulation as in Eq. ( 5.23| ) gives 



21 + 1 100^ ' 



G77oZcqln ( — 



\] n(^)(fc) r 




'J km I 





(6.23) 



(6.24) 



Substituting Eq. (|6.24| ) into Eq. (|6.lD and using Eq. ( p.22| ), we obtain 

2 



EE(T) 



49 
40000 

X 



(27r)2"+i2; 
dk 







1 



Grj^Zcq In 



3 V A:_D 



^cq, 
2n+3 



ikoVo 



i2n+3 



r2(2±3)(2?l + 3) (fcA%)2"+6 
1 



Ji+Ak-qo) + 



2{kr,o) 



kr]o 



(6.25) 



Note that (27r)^"+^"i3|/ [F^ (^^) ^^"^^] oc A '^, wh ere A is the normalization of the magnetic power spectrum in 
Eq. (2^). For n > —3/2, using Eqs. ( ^.8| ) and ( |3.13 ) and keeping only the highest-order terms in I gives 



49 



20000' 



(2tt 



,2n+ll 



G77oZcqln ( — 



[koria 



i2n+3 



r2(«±3)(2n + 3) (fcA77o)2"+e 



In 



kom 



n > -3/2. 
(6.26) 



For —3 < n < —3/2, using Eqs. (5.13) and (5.14) and keeping only the highest-order terms in /, we obtain 

2 



j2^EE{T) _ r,2n-7 ^9 



625^ ' 



G?7o2;oqln ( — 



F(-2n- 3) 



BliAn^ + i) 



r2(-n - 1)F2 (2±3) {2n + 3)(n + l)(n + 3) V^A^yo 



2ji+6 



-3 < n < -3/2. 



(6.27) 



From the properties of radial functions, Hu and White place upper bounds on how fast various power spectra 
can grow with I (see Eq. (78) of [^). In particular, tensor polarization power spectra can grow no faster than 
p^EE,BB(T) ^ results for the tensor E- and B-type (Sec. VI B 2) polarization power spectra seem to violate 

this constr aint f or n > — 2 by an additional factor of /, which arises fr om n umeri cal ap proximations as in the second 
line of Eq. ( |5.23| ). Within the tensor integral solutions of Eqs. ( 5.22| ), ( |6.23 ), and ( 6.3Q| ), we have to evaluate integrals 
of th e for m J^" dx [j2{x)/x][ji{xo — x)/{xq — xY\. The piece j2{x)/x comes from the gravitational wave solution h of 
Eq. ( 4.2l| ) whereas the piece j;(a;o — a;)/ (xq — x)^ comes from the radial fimctions. In Ref. |59|, only the radial function 
properties are used to determine the upper bounds on the power spect ra growth rate, whereas the source behavior 
has been entirely neglected. Our numerical approximation in Eq. ( 5.23| ) takes into account the source behavior, i.e. 
j2{x)/x, and this introduces an additional factor of I in the resulting power spectra. Note that besides the tensor 
polarization power spectra, all the remaining power spectra conform to the growth constraints given by Ref. [p9| . 



2. B-type Polarization 
The B-type polarization integral solution for tensor perturbations is [p9| 
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^ Jo 



21 



where 



3[{x) + 2 



(6.28) 



(6.29) 



is the tensor B-type polarization radial function. Using Eqs. ( p. 20 ), (4.21) and ( |6.5D , and defining x = krj and 
Xq = krjQ, we approximate the tensor B-type polarization integral solution as 



21 + 1 



(27r) 



Gr^o^cq In 



cq 



A similar manipulation as in Eq. ( 5.23 ) gives 



21 + 1 -lOO^^'^)^^ 



n(^)(fc) 



{1 + 2} ji+i[xo-x) 



Gr/gZoqln — 



km 



Xq — X 



i'-i4^-Mkm) 



Substituting Eq. ( 6.31 ) into (6^) and using Eq. ( 2.22 ), we obtain 



BB{T) 
I 



49 
40000 

X 



(2^) 



dk 



2n+12i 







Gr/Q Zcq In ( 



fc 



n + 3 Kko 



^cq 
2n+3 



Bim 



(6.30) 



(6.31) 



(kom) 



2n+3 



r2(ii±3)(2n + 3) (fcA??o)2"+6 

,^i+3(fc??o) 



Ji+i{krio) 



(6.32) 



Note tha t (27r)^"+^°_B^/ [F^ (^^^) ^A^llf ] A ^, wh ere A is the normalization of the magnetic power spectrum in 
Eq. (U). For n > -3/2, using Eqs. and (^.13D , we obtain 



l2(jBB(T) 



49 



;(2^) 



2n+ll 



Gl]lzcq In 



20000 ^ 

n > -3/2. 

For —3 < 71 < —3/2, a similar calculation gives 



cq 



{koTlo) 



2n+3 



r2(«±3)(2„ + 3) (fc;,r7o)2"+6 



H I ) 


- 1 







(6.33) 



p^BBiT) ^ 22«-4i^(27r)2"+12 



625 

-3 < ?i< -3/2 



Gr/gZcqln ( 5^ 



F(-2n - 3) 



-Bin 



F2(-n- 1)F2 (n±3) (2n + 3)(n + l)(n + 3) V^A^yo 



2n+6 



(6.34) 



VII. CROSS-CORRELATION POWER SPECTRA 



Since temperature O; has electric parity (—1)', only Ei couples to 8; in the Thomson scattering and hence Gf ^ is 
the only possible cross correlation. The cross-correlation power spectrum is defined similarly as the temperature and 
polarization power spectra: 

G^-) = 1 (7.1) 
' TT J 21 + 1 21+1 ' ^ ' 

where X stands for V or T. 
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A. Vector Cross-Correlation Power Spectra 



As discussed in Ref. ||5^] and shown in its Fig. 5, the vector dipole radial function j^^^^ does not correlate well 
with its E-type polarization radial function e[^' whereas its quadrupole radial function jp^'' does. Therefore to 

compute the vector cross-correlation power spectra, we need to retain the term proportional to iP^'' in the vector 
temperature integral solution, though in the calculation of the temperature power spectra, we have neglected it since 
it is suppressed relative to the term 



Beginning with Eq. ( p.q ), retaining the jl 



(2V) 



we arrive at the following vector temperature integral solution as in Eq, 



term, neglecting the vector potential, and using Eqs. (^) and (|£ 

6|): 



er^(r7o,fc) 



27 + 1 



7 dec 



(7-1) 



j/(fc'7o) ji+iikfjo) 



krio 3 [ (fc?7o)^ krio 

Substituting Eqs. (7.2) and (6.6) into (7.1) and using Eqs. (4.14) and ( p.l^ ), we obtain 

,,4 



(7.2) 



BE{V) 




(kDVo] 



2n+3 



r2(li±3) (2n + 3) (fcAr7o)2"+6 

Ji+i/2{kvo)Ji+3/2{kvo) 



{kvor 

Ji+i/2{kvo)Ji+3/2{kVQ) 



{km? 



(kvo 



{krio) 



(7.3) 



where again we have neglected the contribution coming from the damped vorticity term, which is negligible (< 3%) 



for 7 < 500 and all cases of n. Note that (27r)^ 



-•AX 



I r (^) k 



2n+61 
A 



cx , where A is the normalization of 



the magnetic power spectrum in Eq. (2_^). The first two terms within the curly bracket arise from correlating jl 



with 



Although the second term proportional to Ji+i/2{k'i]o)Ji+3/2{kriQ) cannot be approximated analytically, 
a numerical evaluation however shows that these two terms always roughly cancel each other, which agrees with 



Ref. 1^ 



that jp^^ does not correlate well with e'l'' ' . The remaining three terms arise from correlating ji^^ ' with 
In the limit 7^1, these terms and the three Bessel terms within the vector E-type polarization power spectrum 



(2V) 



expression of Eq. (6.7) are almost identical. To simplify the approximation, we will neglect the two terms arising from 
the correlation between and e[^^ Thus apart from an overall minus sign, the resulting power spectra fo r al l 



cases here are appr oximately equal to the corresponding vector E-type polarization power spectra, given in Eqs. (6.9) 
and (|6Tl| ) to (|6T3| ). 

Since we have neglected two terms in the vector cross-correlation power spectrum expression of Eq. (^^), accuracy 
here is worse than that of the corresponding E-type polarization power spectra. Note that the terms arising from the 
correlation between jp^"* and e^^"* are suppressed by an additional factor of kL^ dec relative to the two terms arising 
from the correlation between and e|^^. Because of this suppression factor, a numerical calculation shows that 

the residuals of the first two neglected terms can easily be the same order as the remaining retained terms, reducing 
the accuracy of our approximation. Our approximation is good to within a factor of three in general and tends to 
underestimate. 



B. Tensor Cross-Correlation Power Spectra 

Using Eqs. ( p5| ), ( |6^ ), and ( |7l| ), we obtain the tensor cross-correlation power spectrum expression 



CI 



eE{T) 



49 
20000 

X 



(27r 



,2n+12^3 



Grilzcq In 



{kom) 



2ri+3 



r2 (li±3) (2n + 3)770 (fcA?7o)2"+^ 



dk k 







'i\k 



2n+3 



2(A:77o 



Jf+sikVo) 



Ji+3ikVQ)Ji+4{kriQ) 
kT]Q 



(7.4) 



Note that (27r)2"+i°i3^/ [r^ (^^) fc ^"^^J cx A^, where A is the normalization of the magnetic power spectrum in 
Eq. (2.5). For n > —3/2, using Eq. (5.13) and keeping only the highest-order terms in 7, we obtain 



20 



49 
15000 



G?7oZoqln — 



For —3 < n < —3/2, a similar calculation gives 



' 625^ 

-3 < n< -3/2. 



G77oZcq In 



cq 



2n+3 



r2(n±3)(2n + 3) (fcAr?o)2"+6' 



> -3/2. 



(7.5) 



r(-2n-l) B|(2n-1) 
r2(-n)r2 (nfiy (2n + 3)(n + 3) V^aJ^o 



2n+6 



(7.6) 



VIII. RESULTS AND DISCUSSION 



The CMB power spectra generated by a stochastic magnetic field are plotted for ; = 5 to Z = 500 in Figs. 0, |, and|. 
Since we are interested in the signatures of the various microwave background power spectra arising from primordial 
fields that are large enough to result in the observed galactic fields via adiabatic compression, for each plot, we choose 
the magnetic comoving mean-field amplitude to be B\ = 10~^ G and fix A = 1 Mpc, i.e. galaxy and cluster scales. 
For simplicity, we consider a standard Cold Dark Matter Universe (sCDM), i.e. a flat Universe composed of only 
dust and radiation (ryo — 6000/i~^ Mpc) with fli, = 0.05 and h = 0.5. Including a possible cosmological constant will 
affect the scale factor evolution only relatively recently at redshift of a few and will result in a slightly larg er r/o- The 
magnetic power spectrum cutoff wavenumbers for vector and tensor perturbations are given by Eqs. (3.3) and (3.6) 
respectively. Thus for n = —1 and n = 2 for example, with B\ = 10~^ G and A = 1 Mpc, we have ko — 27.9 Mpc~^ and 
ki) ~ 14.7 Mpc~^ respectively for vector perturbations; whereas for tensor perturbations, we obtain k^ ~ 80.9 Mpc~^ 
and ku ~ 27.1 Mpc"-^ respectively. For the tensor perturbations, we assume Zin/^cq = 10^ as in Ref. |5|]; the resulting 
fluctuations, however, depend only logarithmically on z-m- In our analysis, we do not decompose the magnetic fleld 
into a large homogeneous component and a small fluctuating piece. The stochastic magnetic field then affects the 
stress-energy tensor and hence the metric perturbations quadratically. In computing the source terms of vector and 
tensor perturbations [cf. Eqs. (2.18) and (2.22) respectively], convolution of the magnetic field couples the large and 
small scale modes, resulting in the cutoff scale perturbations completely dominate the large scale modes for n > —3/2. 
Thus for n > —3/2, kjj will determine the overall amplitude of the fluctuations. 

Throughout the paper, we have been stating explicitly the terms that are proportional to the normalization A of 
the magnetic power spectrum. Any power-law magnetic field can be specified completely by the normalization A and 
the spectral index n. Since we are interested in constraining the prim ordial magnetic field strength on galaxy scales, 
we choose to fix B\ and A and determine A for each n using Eq. ( p.7| ) . If how ever one is interested in the CMB power 
spectra with A fixed, then each n will give a different value of B\{n) via Eq. (2.7). Either way will not affect the final 
constraints for the magnetic comoving mean-field amplitude. Indeed, we find it easier to constrain the amplitude by 
keeping Bx fixed. 

Figure |l| shows the separate vector and tensor contributions to the CMB power spectra for four different values 
of n. For n < —3/2, the CMB power spectra do not depend on ko, and the size of the anisotropies increases as n 
gets smaller; a scale-invariant magnetic fleld with n = —2.99 generates the largest anisotropies and hence it will yield 
the most stringent limit on the primordial magnetic field (see also Fig. 1 of |5^). For n > —3/2, the CMB power 
spectra are fc^i-dependent and scale as k^^^; thus more and more stringent magnetic field limits can be obtained as 
n increases toward causal value^. For the vector perturbations, the BB power spectrum is slightly larger than that 
of the EE's, as also pointed out in Ref. that the vector CMB polarization is dominated by the B-typc modes; 
whereas the EE and 8E power spectra are approximately identical. Naively, the 8E cross correlation would be 
expected to be larger than the polarization power spectra simply because the temperature fluctuations are larger than 
the polarization fluctuations. However, the temperature fluctuations are dominated by the vector dipole term, which 
correlates poorly with the radial function describing E-type polarization. Thus the 0E spectrum is dominated by a 
subdominant temperature contribution arising from the vector quadrupole term, which then coincidentally renders 
the spectrum a form approximately identical to the E-type polarization itself. However in reality, the 0E spectrum 
can be slightly larger than the polarization spectra since our approximation is good to within a factor of three and 



^Figure 1 of Ref. shows a weaker and weaker upper bound for Bx as n increases from —3/2. This is because the authors 
there have inappropriately adopted a smoo thing scale smaller than the magnetic damping scale and have employed a different 
value for the damping scale (see Sec. IIIB for details). 
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500 



FIG. 1. The microwave background power spectra for vector (left panels) and tensor (right panels) perturbations from a 
power-law stochastic magnetic field with spectral index n. Solid line represents &&, dash-dot line represents EE, dotted 
line represents BB, and dash-dot-dot-dot line represents OE. The magnetic comoving mean-field amplitude is chosen to be 
Bx = 10~^ G, with a smoothing Gaussian sphere comov ing radius of A — IMpc. The magnetic damping cutoff wavenumbers 
for vector and tensor perturbations are given by Eqs. (3.3) and (3.6) respectively. The absolute values of the vector cross 
correlations are plotted. For the tensor perturbations, we assume Zin/zcq = 10^. 
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FIG. 2. Same as in Fig. except that the microwave background power spectra are for vector plus tensor perturbations. 
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FIG. 3. 
represents 
n = 2.0. 



Each panel shows a single power spectrum for various values of n. Solid line represents n = —2.99, dashed line 
n = —2.0, dash-dot line represents n — 0.0, dash-dot-dot-dot line represents n — 1.0, and dotted line represents 
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tends to underestimate in general (see Sec. VII A). Note that while n —3 corresponds to a scale -inva riant magnetic 



field, the vector power spectrum is not flat for this value. The reason is that the vorticity, Eq. ( 4.14 ), has an extra 
factor of k compared to the magnetic field itself. The vector cross correlation is always negative; its absolute value is 
plotted. 

For the tensor perturbations, the E-type is slightly larger than the B-type polarization power spectrum. The 
polarization power spectra are actually comparable to the temperature power spectrum for n > —3/2. This is due 
to the additional l ogari thmic depe ndence on the magnetic damping cutoff wavenumbcr for the polarization power 



spectra [cf. Eqs. (3.26) and (|6.33|)], and also because both the temperature and polarization fluctuations are due 



to the intrinsic temperature quadrupole moments, which arise from the gravitational wave solution h of Eq. (4.21) 
instead of being generated via free streaming the dipolcs as in the case of the vector perturbations. Also for n > —3/2, 
the gravitational wave source term is approximately independent of k and the resulting power spectra then possess 
the well-known behavior of a white noise source I'^Ci oc l^. Furthermore, since tensor perturbations are damped on 



scales smaller than that of the vector perturbations as discussed in Sec. [II B, their induced anisotropies will then be 
larger than that of the vector's for n > —3/2 where ko determines the overall amplitude of the microwave background 
power spectra. As expected, the tensor power spectrum is flat for n — > —3 since we have a scale-invariant magnetic 
field for this value. The tensor cross correlation is always positive. 

The difference between the sign of the vector and tensor cross correlations can be understood from the geometric 
properties of the projection of their corresponding temperature and polarization sources as anisotropies on the sky 
p9[ . The sign of the vector and tensor cross correlations is determined by respectively (cf. Eq. (80) of Ref. [^) 

sgniCf"^^^] = -sgn[p(^)(fF(^))], (8.1) 
sgn[cf ^(^'] = sgn[p(^) (fF(^) - h)] . (8.2) 

The sign of the vector cross correlation is therefore always opposite to that of the tensor cross correlation. 

Each panel in Fig. ^ shows the total vector plus tensor contributions for the various power spectra for a particular 
value of n. Each panel in Fig. || replots one of the four power spectra for a range of spectral indices n. As the 
spectral index becomes greater than zero, the amplitudes become quite large. Again, this is because for n > —3/2, 
the magnetic cutoff wavenumber ko determines the overall amplitude of the power spectra. For a scale-invariant 
n —> —3 spectral index, a magnetic comoving mean-field amplitude of B\ = 10^^ G, and a comoving smoothing scale 
of A = 1 Mpc, at I = 500 for example, the temperature and polarization power spectra are smaller than the amplitudes 
as expected from scale-invariant density perturbations normalized to COBE (i.e., CMB fluctuations in "standard" 
cosmological models). But for n = 0, the temperature power spectrum is essentially a factor of 50 whereas the E-type 
polarization power spectrum is a factor of 10* larger than that expected from the scalar sCDM model at I = 500. 
Therefore, observational limits will be much stronger for causal fields than for scale-invariant fields. 

To estimate the potential observational limits on stochastic magnetic fields, the induced microwave background 
anisotropies must be large enough to be disentangled from the anisotropies arising from density perturbations. Current 
temperature maps give power spectrum measurements with error bars on the order of 10% out to I = 400 for bins 
of width = 50 ||l^,|IJ]. The MAP satellite, which is already in orbit, will make temperature measurements out to 
around I = 800 and will reach the cosmic variance limit, AC; = (^ -I- l/2)~^/^C;, out to I — 400 By the end 

of the decade, and perhaps within five years, we can expect a cosmic-variance limited temperature power spectrum 
measurement to I ~ 3000. Polarization fluctuations will also be detected soon, and the progress in their measurement 
will likely lag temperature fluctuations by about a decade. A rough but conservative estimate is that a magnetic-field 
signal which is at least 10% of the dominant density-perturbation signal will be detectable. The ultimate sensitivity 
in measuring the temperature power spectrum at, say, I = 500 will be significantly better than this, and the extent to 
which magnetic fields can be constrained depends more on the degeneracy of the magnetic field signal with shifts in 
various cosmological parameters. Basic statistical techniques for pursuing such an analysis are well-known (see, e.g., 
p2[) and will be considered elsewhere. 

Using this crude 10% criterion, we can anticipate constraints on stochastic magnetic fields from upcoming temper- 
ature measurements (e.g., the MAP satellite, currently taking data) by simply comparing the predicted amplitude 
at I — 500 to the amplitude of current measurements, which is on the order of PCi ~ 10^^. We assume that the 
remainder of the power spectrum is used for discrimination between the signals from magnetic fields and other tem- 
perature power spectrum contributors. For the scale- invariant magnetic field with n — > — 3, a comoving mean-field 
amplitude of B\ — 10^^ G gives temperature anisotropies at the level of ~ 3% of current measurements. Since for 
n < —3/2, PCi cx B^, the constraint from temperature perturbations on a comoving IMpc scale will be around 
1.4 X 10^^ G, which is approximately at the same level as the previous constraints for a primordial homogeneous 
magnetic field Q-^^. The addition of E-type polarization measurements here will improve the constraint, since 
the ratio of the E-type polarization to temperature power spectra is larger for stochastic magnetic fields than the 
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dominant density perturbations. For n > —3/2, the polarization power spectra are comparable to the temperature 
power spectrum due to the dominant tensor perturbations; thus E-type polarization measurements will yield more 
stringent constraints than temperature measurements alone. Here we will be conservative and project stochastic 
magnetic field constraints using temperature measurements only. For n > —3/2, wc have PC; cx where 
14/(n + 5) = 4 +J-2/(n + 5)](2n + 3), since I'^Ci cx A^fc^^+s, A Bl [cf. Eq. (^j])], and ko oc ^^^/(n+s) j^^^ 
Eqs. ( |3.3|) and (3.6)]. As n increases towards causal values, the amplitude of the temperature fluctuations increases 
as k^j^^'^ and hence the constraints become stronger. At n = 0, PCi at I = 500 is approximately 5 x 10^^, which will 

yield a constraint on Bx of I q^^^^ "q-'j ) x 10"^ G ~ 10^^° G. For the causal field n — 2, the constraint on B\ 

will be as small as 4 x 10"^'^ G. Such constraints will be stronger than any current limits on Mpc-scalc primordial 
stochastic magnetic fields at decoupling. 

Ultimately, B-type polarization has the greatest potential for constraining primordial magnetic fields. This is 
a cleaner signature, because primordial scalar (density) perturbations produce none p0||7^ . Aside from polarized 
foreground emission, the only other expected sources are from primordial tensor perturbations and from gravitational 
lensing |Q . Tensor perturbations with a spectrum near scale-invariant will give significant anisotropics only at large 
angular scales [l < 100), while lensing contributes mainly at small angular scales (/ > 500). Stochastic magnetic fields 
will contribute on intermediate scales and should be clearly distinguishable. If foreground emission can be separated 
from its frequency dependence, limits on B\ from B-type polarization should be determined purely by measurement 
error bars on Cf^ . Note that a primordial magnetic field also generates an additional B-type polarization signal 
via Faraday rotation of the CMB polarization |^^. This signal will be negligible compared to the direct B-type 
polarization signal for any frequency of practical interest. 

All of the results in this paper have been obtained via analytic approximations to the exact solutions. Apart from 
the vector cross-correlation power spectrum, the accuracy of the results is as good as the quality of the analytic 
approximations to various expressions, except that the vector temperature case has neglected an additional few 
percent temperature contribution arising from the angular dependence of polarization and the vector potential. These 
approximations are all discussed in the text; in sum, they are good to within 20% over the range of parameters 
considered, with the exception of the vector temperature case in the regime —2.5 < n < —2, which is good to within 
30%. Meanwhile, accuracy of the vector cross -cor relation power spectrum is only good to within a factor of three, 
since we have neglected the two terms in Eq. (7.3) arising from the correlation between the temperature dipole and 
the E-type polarization radial functions. It is important to realize that errors in these analytic approximations will 
have negligible effects on the estimation of the magnetic field limits for n < 2 since the amplitude of each power 
spectrum scales as Bj^ for n < —3/2 and for n > —3/2. 

In this paper, we have focussed on the magnetic field-induced microwave background anisotropics for I < 500, 
where the analysis is relatively simple and free from the detailed microphysics of recombination. We have only 
considered vector and tensor metric perturbations; for smaller angular scales 500 < I < 2000, the magnetic- induced 
CMB anisotropics are dominated by vector perturbations [^. Stochastic magnetic fields will also produce scalar 
perturbations. This case is significantly more complex due to physical compensation effects and the large number of 
terms involved in the relevant expressions. Rough analytic estimates show that including the scalar results will only 
modestly improve the magnetic field constraints given in this paper, since radiation pressure prevents the induced 
density fluctuations from growing effectively before recombination and the compressional modes are erased up to the 
Silk scale Ls |41 4^. Therefore scalar perturbations will generally give a subdominant contribution to the microwave 
background anisotropy. 

Our results suggest that while it may be plausible for primordial stochastic fields with n ^ to result in the observed 
galactic fields via adiabatic compression alone, it will be very difficult for causal fields without invoking some form 
of dynamo mechanism. In a recent paper | |7^ , a similar calculation using the nucleosynthesis bound on gravitational 
radiation induced by the anisotropic stress of a primordial stochastic magnetic field yields extremely stringent limits 
on the galactic-scale magnetic field amplitudes {B\ < 10~^^ G) for fields generated at the electroweak phase transition 
or earlier, thus ruling out most of the magnetogenesis processes for primordial fields seeding the observed large-scale 
coherent galactic fields. 
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APPENDIX A: DERIVATION OF THE VECTOR ISOTROPIC SPECTRUM 

Our objective is to derive the vector isotropic spectrum |n'^^^(fc)p defined in Eq. (2.17), which will be useful for 
calculating vector CMB power spectra. Using Eq. ( p. 14 ), the two-point correlation function of n^^-* is given by 

^n(^)(k)n(^)*(k')) - p,,fc./^,fc:(r(f (k)r(f *(k')), (Al) 



where P'^^ — 5id — k[k'j^. We simplify our calculation by splitting the electromagnetic stress-energy tensor into two 
pieces: T^P^k) = r^^'^^k) +T^^'^\k) where 

1 1 



^ '~ (27r)3 4tt 

^. ^(k) = (5. 



(27r)3 Stt 



J d^pB,ip)B,ik-p), 

J d^pBiip)Bi{k-p). 



(A2a) 
(A2b) 



The two-point correlation function of the electromagnetic stress-energy tensor in Eq. (Al) will now be described by 
a sum of four two-point correlation functions: 



(k)rif *(k')) = (rif ^^^(k)r(-^)*(k')) + (rff ^^)(k)r(f-)*(k')> 

+ (-r^(k)rif^)*(k')) + (r(f^^(k)r(f^)*(k')). 



(A3) 



Lib ^cd 



above has a non-vanishing contribution toward the two-point correlation function of H^^^ in 



Only 

Eq. (Al), since each of the remaining correlation functions in Eq. ( |A3D contains cither Sat, 5cd, or both, and will 
vanish when they are acted upon by PujkaP^Jt'c- We can now rewrite Eq. (Al) as 



nr^(k)nr^*(k')) = p.bkaPUK{rir'{^)<f'^*{^))- 



(A4) 



We can e valuate the two-point correlation function {t^^ ' t^^ '^'") as follows. Beginning with the definition of 
Eq. ([A2a|), we assume the random magnetic field is Gaussian and apply Wick's theorem 



{B,{k,)B,{kj)Bi{ki)B^{k^)) = (B,{k,)Bj{k,)){Bi{ki)B^{krn)) + (B.(k,)Si(kO)(B,(k,)B„(k™)) 

+ (B,(k,)S™(k™))(B,(k,)Bi(kO) 



(A5) 



and the reality condition B*{k) = Bi{—k), and finally use Eq. ( pH] ) for the form of the stochastic magnetic field 
two-point correlation function to arrive at (see also pm) 



(-lf'^^(k)rif^^k')) = 



1 



(47r)2 



d^pP{p)P{\k - p\)[{Sac - PaPcKSbd - (k - p)fc(k - p)d) 



{Sad - PaPd){Sbc - (k - p)h(k - p)c)]<5(k - k') 



(A6) 



Substitute Eq. ( |A6|) into Eq. ( |A4| ) and define 7 = k • p, /3 = k • (k — p), and ^ = p • (k — p) to obtain the two-point 
correlation fmiction of the vector n,-^-*: 



nr)(k)nr^*(k')) 



1 



(47r)2 



d'pPip)P{\k - p|)[(l - 72)(1 + /32) + jf3(^ „ ^f3)]6{k ~ k'). 



(A7) 



The integral above is similar to the mode-coupling integral P{k) in Eq. (11) of Ref. Although it cannot be 

evaluated analytically, terms within the square bracket are products of cosine factors; hence the s quare bracket itself 
can be approximated by unity, which has essentially been done in Ref. js^. Comparing with Eq. ( 2.17| ) gives 



(A8) 
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Using the expression for P{k) in Sec. || and choosing k to be the polar axis, the vector isotropic spectrum becomes 



|n(^)(fc)| 

The integral over 7 is 



Bj 

32 r2 (Ii±3) fc2„+6 



1/2 



1/2 _ 



1 



[{k+pr+'-\k-p\ 



n+2l 



^1 A:p(n + 2) ' 

and the expression within the square bracket above can be approximated as 



{k + pT+^ -\k-p\ 



„+2_ /2(n + 2)A:"+V, P<k- 

2(n + 2)/cp"+^, otherwise. 



(A9) 



(AlO) 



(All) 



Substituting Eqs. ( |A10[ ) and ( |All|) into Eq. ( [A9| ) and evaluating, we finally arrive at the expression for the vector 
isotropic spectrum, Eq. (^.18). 
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